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Angular power spectrum of gamma-ray sources for 
GLAST: blazars and clusters of galaxies 
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ABSTRACT 

Blazars, a beamed population of active galactic nuclei, radiate high-energy 7-rays, and 
thus are a good target for the Gamma Ray Large Area Space Telescope (GLAST). As 
the blazars trace the large-scale structure of the universe, one may observe spatial 
clustering of blazars. We calculate the angular power spectrum of blazars that would 
be detected by GLAST. We show that we have the best chance of detecting their 
clustering at large angular scales, 9 > 10°, where shot noise is less important, and the 
dominant contribution to the correlation comes from relatively low redshift, z < 0.1. 
The GLAST can detect the correlation signal, if the blazars detected by GLAST trace 
the distribution of low-z quasars observed by optical galaxy surveys, which have the 
bias of unity. If the bias of blazars is greater than 1.5, GLAST will detect the corre- 
lation signal unambiguously. We also find that GLAST may detect spatial clustering 
of clusters of galaxies in 7-rays. The shape of the angular power spectrum is different 
for blazars and clusters of galaxies; thus, we can separate these two contributions on 
the basis of the shape of the power spectrum. 

Key words: gamma-rays: theory — BL Lacertae objects: general — galaxies: active 
— galaxies: clusters: general — cosmology: theory — large-scale structure of Universe. 



1 INTRODUCTION 

Active galactic nuclei (AGNs) are highly energetic astro- 
physical objects, which are often accompanied with rela- 
tivistic jets powered by accretion onto supermassive black 
holes located at the central region of galaxies. The AGNs 
radiate in a wide frequency range, from radio waves to 
7-rays. The Energetic Gamma Ray Experiment Telescope 
(EGRET) on board the Compton Gamma Ray Observatory 
(CGRO), which detects 7-rays in GeV energy, have found 
~ 60 AGNs with its all-sky survey campaign, and all of them 
but o ne (M87) were classified as "blazars" (IHa rtman et al.l 
Il999h . The features in the spectrum and light curves of these 
blazars indicate that the relativistic jets are directed towards 
us. 

The Gamma Ray Large Area Space Telescope (GLAST) 
is equipped with the Large Area Telescope (LAT) instru- 
ment, which is an upgraded version of the EGRET. Its 
large effective area (10 4 cm 2 ) as well as good angular reso- 
lution improve a point-source sensitivity by almost two or- 
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ders of magnitude compared to the EGRET, which would 
increase the source statistics significantly. It is expected 
that a thousand to ten thousand bl azars would be de- 
tecte d as point sources _by GLA S T Istecker fc Salamonl 
19961: IChiang fc Mukherjeel Il998l : iMiicke fc Pohll |2000| ; 



Narumoto fc Totanil l2006t ). Such a dramatic improvement 



in point-source sensitivity would allow us to determine the 
7-ray luminosity function (GLF) of blazars with unprece- 
dented accuracy. In contrast, the current blazar GLF is 
based upon merely ~ 60 blazars detected by EGRET. 

As the blazars should trace the large-scale structure 
of the universe, they should exhibit spatial clustering. In 
this paper, we investigate whether the spatial clustering of 
blazars is detectable by GLAST, especially by focusing on 
the angular power spectrum, a quantity projected along the 
line of sight. This would provide the first direct measurement 
of the bias parameter of blazars, immediately after GLAST 
starts taking the data, which should be compared with the 
results from other classes of AGNs in order to get further 
insight into the unification picture of AGNs. In addition to 
blazars, it has been pointed out that GLAST may detect 7- 
ray emission from clusters of galaxies l|Colafrancesco fc Blasil 
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ll99Sl ; lTotani fc KitavamalbOOOt ). We thus also calculate the 
spatial clustering of clusters of galaxies in 7-rays and discuss 
its detectability with GLAST. 

This paper is organized as follows. In Section [2] we cal- 
culate the angular power spectrum of blazars that would 
be detected by GLAST. In Section[3]we study detectability 
of the spatial correlation of blazars and discuss current (in- 
direct) observational constraints on the bias of blazars. In 
Section|4]we calculate the angular power spectrum of galaxy 
clusters. Section [5] is devoted to further discussions, and we 
conclude in Section [6] 



2 ANGULAR POWER SPECTRUM OF 
BLAZARS 

2.1 Formalism 

The angular power spectrum of blazars that would be de- 
tected by GLAST is given by the sum of the shot ( Poisson) 



noise term, Cf, and the correlation term, Cf , as l|Peeblesl 



1980) 
Ci 

cf 



cf + cf, 



AT 1 , 



dcosfl Pi (cos B)w (8), 



(1) 
(2) 

(3) 



where J\f = dN/dQ is the number of blazars per solid angle, 
and w(9) is the angular correlation function of blazars that 
would be detected by GLAST. Note that the shot noise term 
is independent of multipoles. 

A standard procedure to calculate the angular correla- 
tion function is as follows. We model the 3-d spatial correla- 
tion function of blazars, £(r, z), as the correlation function of 
dark matter particles, multiplied by the "bias" factors that 
depend on physics of formation and evolution of blazars in 
dark matter haloes. We then project the resulting 3-d cor- 
relation function on the sky to calculate the 2-d angular 
correlation function of blazars, w(9). As the bias factors de- 
pend on redshift and luminosity of blazars, we model £(r, z) 
as £(r; L 7 ,i, L 7 , 2 |z) = £un(r, z)&g(L 7 ,i, z)& s (L 7j2 , z), where 
r — \x2 — x\\ is the distance between two blazars, L 7j i and 
L 7i 2 are their luminosities, and fii a (r, z) is the 3-d correla- 
tion function of linear dark matter fluctuations. As we show 
in this paper the angular correlation function of blazars may 
be detectable only on large scales, and thus the linear cor- 
relation function and the linear bias model would provide a 
good approximation. By proj ecting the 3-d correlation func- 
tion on the sky, one obtains (|PeeblesHl980h 



N 2 w{6) 



dz 



d 2 V 
dzdQ 



X(z) <t>(z) 2 b B (z) 



X J dutiun(Vu 2 +X(z) 2 2 ,z) 



(4) 



where x( z ) i s the comoving distance out to an object at z, 
d 2 V/dzdfl is the comoving volume element per unit solid 
angle and per unit redshift range, bs(z) is the average bias 



of blazars weighted by the 7-ray luminosity function (GLF) 
of blazarsQ p 1 (L 1 , z): 



&s ( 2 ) = TTt/ di 7 p 7 (L 7 ,z)6 B (L 7 ,z), (5) 

<P( Z ) iL 7 (F Tilim ,z) 

and 4>{z) is the cumulative GLF of blazars, i.e., GLF inte- 
grated from a given minimum luminosity cut-off, 



dL 1 p 7 (L 7 , z). 



(6) 



-M-FVlim.z) 



Note that we have not used the so-called "small-angle 
approximation" or "Limber's approximation," as we are 
mainly interested in the signals on large angular scales, 
9 > 10°. 

We calculate £n n {r, z) from the power spectrum of linear 
matter density fluctuations, Pu n (k): 



6m(r,z) = 



k 2 dk „ ,, , sin kr 
fiin(fc)- 



(7) 



2-7T 2 kr 

We us e the linear transfer function given in lEisenstein fc Hul 
|l999l ) to compute Pi in (/c). 

Equations ([2]) and ((3]) suggest that Cf = A/ --1 dom- 
inates when the number of blazars detected by GLAST is 
small, making it difficult to detect the correlation term. It is 
therefore very important to understand how many blazars 
one can detect with GLAST. In the next subsection we cal- 
culate the expected number co unt of blazars for GLAST 
using the latest GLF of blazars l|Narumoto fc To tani 2006). 

2.2 Gamma-ray luminosity function of blazars 

The basi c idea behind the model of t he GLF of blazars pro- 
posed bv lNarumoto fc Totanil (|2006t ) is that the jet activity 
that powers 7-ray emission from blazars must be related 
to accretion onto the central black holes, from which X-ray 
emission emerges; thus, the X-ray and 7-ray luminosity of 
blazars must be correlated. We use the following relation 
between GLF of blazars, p-,, and X-ray luminosity function 
(XLF) of AGNs, p x - 



p 7 (L 7 



IT \ 

K—p x (Lx,z) 



(8) 



The advantage of this method is that the XLF has been 
determined a ccurately by the extensive study of the X-ray 
background (|Ueda et al.l 120031 ; lHasinger et all 120051 ). and 
thus the predicted GLF would also be fairly accurate, pro- 
vided that the 7-ray luminosity and X-ray luminosity of 
blazars are tightly correlated. Since not all AGNs detected 
in X-rays are blazars, we have introduced a normalization 
factor, «. We relate the 7-ray luminosity, L 7 , and X-ray 
luminosity, Lx , of blazars by a linear relation with the con- 
stant of proportionality given by 10 q : 



W q L } 



(9) 



where L 7 represents vL v at 100 MeV, and Lx is the X- 
ray luminosity integrated over the R OS A T band, 0.5-2 keV. 
(Both are evaluated at the source rest frame.) We convert 
the measured flux to the rest-frame luminosity by specifying 



1 The luminosity function represents the number of sources per 
unit comoving volume and unit luminosity range. 
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or 7 

of ct 7 = 2.2 (Srccku mar et alii 19981 ), while for X-ray ax = 2 



the spectru m of sources: for 7-ray we use an spectral index 
of a-y = 2.2 (Sreekuma 
ijHasineer et al.ll2005h . 

The AGN XLF, px, is given by a doub le power-law in 
lumin osity with an evolution factor f(Lx , z) (Ha singer et al.l 
120051 ): 



px(L x ,z) 

where 
f(Lx,z) 



Axf{L x ,z) 
(In W)L X 



Lx 
L* x 



+ 



L* x 



(10) 



{l + z) p i 



< Zc(Lx)], 



r ~\ P2 

f[L x ,z c (Lx)} [ 1+2 1 +l x) j " [z > z c {L x )l 

where z c is the redshift of evolutionary peak given by 



z c (Lx) 



(L x > L a ), 
{Lx < L a ), 



(11) 



(12) 



and pi and P2 are given by 

Pi = P*i + Pi [log(Lx/erg s" 1 ) - 44] , (13) 
P2 = P*2 + P2 [log(L x /erg s" 1 ) - 44] . (14) 

lHasinger et~ai1 (|2005h have found A x = 6.69 x 1CT 7 
Mpc" 3 , log(L^/erg s" 1 ) = 43.94 ± 0.11, z* = 1.96 ± 0.15, 
log(L a /erg s" 1 ) = 44.67, a = 0.21 ± 0.04, p\ = 4.7 ± 0.3, 
p* 2 = -1.5 ± 0.7, /3i = 0.7 ± 0.3, = 0.6 ± 0.8, 72 = 
2.57 ± 0.16, and 71 = 0.87 ± 0.10. We call this model the 
"luminosity-dependent density evolution" model, LDDE. 

How robust are our predictions from this model? The 
most important parameter for our purpose in this paper is 
the slope of XLF in the faint end, 71, as the expected num- 
ber count of blazars that would be detected by GLAST is 
sensitive to how ma ny blazars there are i n the faint end of 
luminosity function. iNarumoto fc Totanil |2006l) have fitted 
the GLF of blazars detected by EGRET in order to find 
71, q — log(L 7 /Lx) and k, with the other parameters fixed 
at the best-fitting values from the XLF given above. (The 
blazar sample from EGRET was constructed by requiring 
that EGRET sources were identified as blazars by radio ob- 
servations. The probability that the blazars giving the flux 
above the EGRET point source sensitivity also gives the suf- 
ficient radio flux was taken into account in their analysis.) 
They have found that (-yi,q,n) = (1.19,3.80,5.11 x 10~ 6 ) 
best describes the GLF of EGRET blazars. This 71 is larger 
than that obtained from the XLF, 71 = 0.87±0.10, at the 3- 
level, which may imply that a better model is needed; how- 
ever, we do not investigate this point any further and sim- 
ply accept 71 = 1.19 as the canonical value for the GLF of 
blazars. One should come back to this point, however, when 
GLAST flies and collects many more blazars than available 
now. 

One can also calculate the contribution to the ex- 
tragalactic gamma-ray background (EGRB) from blazars 
once the GLF is speci fied. The EGRB intensity has been 
meas ured by EGRET ijSreekumar et al.lll998l : IStrong et aD 
2004)0 and the best-fit model with 71 = 1.19 accounts 
for 25-50% of the EGRB intensity, depending on the as- 



sumed minimum 7-ray luminosity of blazars, L 1 



= 10* 



Table 1. Parameters of the LDDE GLFs and the expected num- 
ber, N, and the surface number density, J\f, of blazars that would 
be detected by GLAST. We have assumed that no blazars fainter 
than L 7iln i n = 10 41 erg s _1 would exist, and GLAST LAT can 
detect the flux down to F~ t n m = 2 X 10 — 9 cm -2 s _1 for 2 years 
of all-sky observations. 



Model 


(q, 71) 


K 


N 


N (sr- 1 ) 


LDDEl^ 


(3.80, 1.19) 


5.11 x 10~ 6 


3100 


250 


LDDE2^ 


(3.80, 1.31) 


3.90 x 10- 6 


6500 


520 



a Best-fit model of the EGRET blazar distribution. 
^ A model explaining 100% of the EGRB intensity. 



to 10 



43 erg s . Here, we assume that no blazars fainter than 



the minimum luminosity would exist. On the other hand, 
blazars can still account for all the EGRB intensity, if the 
blazars can be as faint as Z/ 7 , m i n = 10 41 erg s" 1 , and the 
faint end of the GLF is slightly steeper, 71 = 1.31, than 
the canonical model. The other parameters are given by 
(q, k) — (3.80, 3.9 x 10 -6 ). This model appears to be a bit ex- 
treme, as 71 = 1.31 is inconsistent with the X-ray determina- 
tion, 71 = 0.87T0.10, at the 4.4-<r level. Nevertheless, we use 
this model to show the uncertainty in our predictions from 
the uncertainty in the faint-end of the GLF. Henceforth we 
call the canonical model (71 = 1.19) the "LDDE1" model, 
and the latter model (71 = 1.31) the "LDDE2" model. For 
both models, we adopt L 7iin i n = 10 41 erg s~ 4 as the lower 
luminosity cutoff. 



2.3 Survey parameters, number count, angular 
correlation function, and power spectrum of 
blazars from GLAST 

The flux sensitivity for point sources of the GLAST LAT 
is -F 7 ,ii m = 2 x 10~ 9 cm -2 s _1 for 2 years of full-sky obser- 
vations and for sources having the energy spectral index of 
2; we adopt this value of the flux sensitivity in the follow- 
ing discussions, unless otherwise stated. The 7-ray flux, F 7 , 
represents the flux integrated above i? m m = 100 MeV, and 
it is related to the 7-ray luminosity through 



L~,{Fj,z) = 



47rd|(ci 7 - 1) 
(1 + z) 2 - a ~> 



E m 1 n F, 



(15) 



where cLl is the luminosity distance. One can calculate the 
number of blazars that would be detected by GLAST from 



N = Q 



dz 



d 2 V 
dzdQ, 



(*), 



(16) 



But these estimates are still controversial llKeshet et al . 2001b). 



where we use z max = 5, Q is the solid angle covered (Q — 
An sr for the all-sky survey), and <j>(z) is the cumulative GLF 
given by equation ©. 

For the canonical GLF model that accounts for 25-50% 
of the EGRB intensity (LDDE1) and the lower luminosity 
cutoff of L 7jm i n = 10 41 erg s _1 (hereafter, we use this value 
unless otherwise stated), we obtain ~ 3100. For the GLF 
model that accounts for all the EGRB intensity (LDDE2) 
and the same luminosity cutoff, we obtain N ~ 6500. These 
results are summarized in Table Q] Therefore, it is easier to 
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Figure 1. Redshift distribution of blazars that would be detected 
by GLAST, for LDDE1 (solid) and LDDE2 (dashed) models. (See 
Table [T] for the model parameters as well as for the expected 
number of blazars.) The thick solid line shows the LDDE1 model 



with 



10 4 



J 7,min x< ~' erg S 

LDDEl model with L. 



the larger the £ 7! min 
detected. 



, while the thin solid lines show the 
and 10 42 erg s — 1 : 
blazars would be 



7 , mm = 10 40 erg s 1 
is, the fewer the low- 2 



detect the spatial clustering of blazars in the LDDE2 model 
than in the LDDEl model. 

Figure Q] shows the redshift distribution of GLAST 
blazars predicted from the LDDEl and LDDE2 model. For 
both cases, the distribution exhibits a sharp cut-off around 
z — 0.01, which is due to our assumption that no blazars 



fainter than L 7 



would exist; a larger I/ 7 , m in results in a 



larger cut-off redshift. Nevertheless, since only a small frac- 
tion of the distribution is eliminated by this effect, the to- 
tal number of blazars that would be detected by GLAST, 
N, hardly changes; for example, we expect 3200 and 2900 
blazars to be observed by GLAST for L 7 , m in = 10 40 and 10 42 
erg s _1 (both for the LDDEl parameters), respectively. On 



the other hand, we shall show in Section T3. II that L 



has 



an important consequence for detectability of the anisotropy 
signal. 

Figure [2] shows the angular correlation function, w{ff) 
(left panels), and the correlation term of the angular power 
spectrum, l(l + l)Cf /2ir (right panels), divided by the aver- 
age bias squared, for the LDDEl (top panels) and LDDE2 
(bottom panels) model. In each panel we vary the GLAST 
LAT point-source flux sensitivity, i^iim, from 2 x 10" 



to 



4 x 10 9 cm 1 s i . As expected, the clustering is stronger 
when more sources are observed, i.e., LDDE2 and lower 



3 DETECTABILITY OF THE BLAZAR 
CORRELATION 

3.1 Signal-to- noise vs blazar bias 

As the correlation function and power spectrum are propor- 
tional to the average bias squared, w{6) oc b B and Cf <x b B , 
whether or not one can detect the angular clustering of 
blazars crucially depends on b B - Before we investigate a 
model of the blazar bias, let us ask this question, "how large 
&s should be, in order for Cf to be detected by GLAST?" 

The statistical error in the measurement of Ci is given 
by the following argument. Assuming statistical isotropy 
of the universe, we have 21 + 1 independent samples of 



| aim | (with different m's) per multipole. Here, 



air. 



is the spherical harmonic coefficient of the distribution of 
blazars on the sky. One may thus estimate Ci from Ci = 
J2m=-l l°im| 2 /(22 + 1). The error in Cj is given by 

2Cf 2(Cf + C?) 2 



{5Cl) (2/ + l)AZ/ sky (2/ + l)AZ/ sky (17) 

where Al is the bin size in I space and / sky is a fraction of 
the sky covered by observations. For the all-sky survey like 
GLAST, we may adopt / sky = 1; we note that the point 
source sensitivity becomes worse near the galactic plane be- 
cause of strong galactic foreground. As Cf = A/" -1 is inde- 
pendent of I and depends only on the inverse of the surface 
density of blazars, one can fit it and subtract it from the 
measured C;, leaving only Cf . The error in Cf, however, 
still contains the contribution from Cf . This shows why it 
is important to detect as many blazars as possible (and thus 
reduce Cf as much as possible), in order to measure Cf. 

Figure [3] shows the l-cr error boxes binned with Al — 
0.51 for the LDDEl and LDDE2 model. We show the errors 
for the average bias of 6s = 1 and 3. (Note that we have 
ignored the redshift dependence of 6s-) We find that it would 
be difficult to detect Cf for the LDDEl plus 6s = 1 model, 
while the other models yield sufficient signal-to-noise ratios. 

To increase statistical power one may sum C; over mul- 
tipoles. Let us define the angular power spectrum averaged 
over 2 «S I 300 



1 30 

C(2^^30) = -^C, 



29 

1=2 

The errors of this quantity is then given by 

30 



(18) 



(scy 



l)f 



1 

29 s " ^ 



(21 + l)/ sky 



(Cf + Cf) 2 . 



(19) 



Figure [4] shows C (2^1^ 30) as a function of the average 
blazar bias, 6s, for the LDDEl (top panel) and LDDE2 (bot- 
tom panel) models. The expected l-cr errors as well as the 
— p 

Poisson contribution, C , are also shown. For the LDDEl 
model we find that GLAST can detect C° if 6s > 1.2. For 
the LDDE2 model the detection is much easier, even for 
6s > 0.5. 



3 A dipole component, C\, depends on Earth's motion and is not 
considered here. 
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Figure 2. (Left) Angular correlation function, w(8), and (Right) correlation term of the angular power spectrum, 1(1 + l)Cf /2ir, of the 

—2 

blazars that would be detected by GLAST. Both have been divided by the average bias squared; thus plotted quantities are w(9)/b B and 
1(1 + 1)C^ /(2-Kb B ). The dotted, solid, and dashed lines show the predictions for the limiting flux of -F 7i u m = 2, 3, and 4 X 10 — 9 cm -2 s _1 , 
respectively. The top panels are for the LDDE1 model, while the bottom panels are for the LDDE2 model. 



Our results depend on the luminosity cutoff of the GLF, 
£7, mm, as the correlation at large separations (I < 30) is 
dominated mainly by relatively nearby (less bright) sources. 
We have therefore performed the same calculations with 
different L 7jm i n (with the other parameters of the LDDE1 
model held fixed), and found that the correlation would be 

Q 

detectable (i.e., C /5C > 1) for the average bias greater 
than 0.9 and 1.7, for 
spectively. 



= 10 4U and 10 42 erg s" 



One may also ask how these results would change, if 
we chose other GLF models. The "pure-luminosity evolu- 
tion" (PLE) model has been used traditionally in the litera - 
ture (|Stecker fc Salamon|[l99r3 ; IChiang fc Mukherieelll998l) . 
while t he LDDE model fits the E GRET blazar properties 
better (|Narumoto fc Totanil l2006t ). Motivated by the cor- 
relation between radio an d 7-ray luminosities of blazars, 
IStecker fc Salamonl (|l996t ) used the PLE model to obtain 
the GLF of blazars. We find that the large-angle correlation 
(I ^ 30) is more difficult to detect in their model: the correla- 
tion would be detectable only when 6s > 4.2. Their model, 
however, was not intended to reproduce the redshift and 
luminosity distributions of the EGR ET blazars, and thus 
their fit to these data is not very good. lchiang fc Mukhcricc 



ljl998l ) improved the PLE model by adjusting a few pa- 
rameters such that the model can reproduce the distribu- 
tion of EGRET blazars. (Although the authors did not use 
the radio and 7-ray lum inosity relation, we incorpo rate this 
in our calculations; see iNarumoto fc Totanil (|2006h for de- 
tails.) Again, we find that the correlation signal is more 
difficult to detect in the best-fit PLE model: the correla- 
tion would be detectable only when ha > 6.9. These results 
are because the PLE model predicts the blazar distribu- 
tion that is much more biased towar d the high redshift (see 
Fig. 11 of lNarumoto fc Totanil l|2006f )L and hence, the large- 
separation power (due mainly to low-redshift blazars) is sup- 
pressed. In fact, the results improve if we instead adopt the 
smaller separation, 30 ^ I ^ 300, where the high-redshift 
contribution becomes larger. The sensitivity to the bias pa- 
ramete r goes down to ha > 2.4 and fcg > 3.0, respectively 
for th e lStecker fc Salamonl l|l996h and lChiang fc Mukherieel 
1 19981 ) models. On the other hand, as we have already shown, 
the latest GLF from the LDDE model, which best describes 
the distribution of EGRET blazars, predicts that the corre- 
lation would be detectable for bs of order unity. 
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1 10 100 1000 1 10 100 1000 



Multipole I 

Figure 3. Angular power spectrum of GLAST blazars, 1(1 + l)C;/27r. The dotted lines show the shot noise term, Cf [equation while 
the thick solid lines show the correlation term, CP [equation for the following models: (a) LDDE1, bs = 1, (b) LDDE1, bs = 3, (c) 
LDDE2, b B = 1, and (d) LDDE2, b B = 3. The boxes show the 1-cr errors in Cf binned with Al = 0.51 [equation tTTH . 



3.2 Modeling blazar bias 

The bias of blazars is not known, and thus GLAST may pro- 
vide the first determination of the blazar bias, if it is greater 
than 1.2 (for LDDE1; 0.5 for LDDE2). In this section we 
estimate the blazar bias from the existing observations us- 
ing three different (indirect) techniques: the bias of quasars, 
X-ray observations, and halo model. However, we emphasize 
that none of these estimates can be very accurate (and in 
fact, they disagree each other) — we will need GLAST to give 
us the answer. 



3.2.1 Bias from optical quasar surveys 

If the blazar is truly a beamed population of AGNs, then 
its bias should be correlated with that of AGNs observed in 
other wavebands. It is therefore natural to use such infor- 
mation to estimate the bias of blazars. 

The optical quasar surveys, such as the Two-degree 
Field Quasar Redshift Survey (2QZ) and the Sloan Digi- 
tal Sky Survey (SDSS), may be the most efficient way for 
doing this. Figure [5fa) shows the spe ctroscopic result that 
the 2QZ suggests 



Croom et al]|2005h 



which is al so consistent with the photometric result from 
the SDSS (|Mvers et all 120061 ). By comparing b Q (z) with 
the bias of dark m atter haloes (which is discussed later), 
ICroom et all ll2005h found that the mass of dark matter 
haloes that host quasars is typically ~4x 10 12 Mq, almost 
independent of redshift. 

In order to make a quantitative comparison between the 
results obtained here and those in the previous section, we 
define an "effective bias," b c B (I), by 



6b (0 



Ci (model bias, &_b(z)) 
Ci(b B = 1) : 



(21) 



b Q {z) = 0.56 + 0.289(1 + ^) 



(20) 



where the numerator is Cf calculated with 6s (z) = bq(z), 
and the denominator is that calculated with T>b = 1. If b c g(l) 
is greater than 1.2 and 0.5 for the LDDE1 and LDDE2 mod- 
els, respectively, GLAST can detect Cf from this popula- 
tion. The top curves of Fig. [6] show b c g(l) as a function of 
I. We find that b c g(l) ~ 0.8 for a relevant range of multi- 
poles; thus, Cf is detectable for the LDDE2 model but not 
for the LDDE1 model. We also find that bf(l) increases as 
I does (haloes at higher redshifts are contributing more to 
the small angular scales), although the dependence is only 
modest and is not able to bring the bias to high enough 
values for detection for the LDDE1 model. 
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Figure 6. The effective bias, V§(1), for the LDDE1 (solid) and 
LDDE2 (dashed) models. The top curves are from bg(z) = bn(z), 
while the bottom curves are from the halo model with Lj—M 
relation [equation (122 H . 
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Figure 5. (a) Quasar bias, 6q(z), from the Two-degree Field 
Quasar Redshift Survey, and (b) contribution to the angular cor- 
relation function per In z, dw(8)/[w(8)dlnz], at 9 = 20°, for the 
LDDE1 (solid) and LDDE2 (dashed) models. We have used b Q (z) 
for the average blazar bias, fog(z). 



Since the multipole I is roughly related to the angular 
separation via 9 ~ 180° /I, the angular power spectrum at 

I ~ 10 contributes most to the angular correlation function, 
w(6), at 9 ~ 20°. In Fig.[SJb), we show that the contribution 
to w(9) per lnz peaks at z ~ 0.01 with a tail extending up 
to z ~ 0.1. At such low redshifts the value of the averaged 
bias is ~ 0.8, almost independent of z. 

3.2.2 Bias from X-ray point-source surveys 

The GLF proposed by iNarumoto fc Totanil (|2006l ) that we 
are using in this paper was derived on the basis of a cor- 
relation between 7-ray and X-ray luminosities of emission 
from blazars. Therefore, the bias derived from the spatial 
clustering of AGNs detected in X-ray surveys may provide 
a useful information regarding the bias of blazars. 

Nevertheless, the clustering of AGNs determined from 
X-ray surveys are somewhat controversial. While both the 
angular and 3-d correlation function of the X-ray bright 
AGNs detected by the ROSAT survey s are consistent with 
those from the optical q uasar surveys (jVikhlinin fc Formanl 
1 19951 ; iMullis et all 12004 ). those from Chandra and XMM- 
Newton suggest that these sources ar e clustered more 
strongly than optical l y selected quasars jYang et al.ll2003l ; 
iBasilakos et al.l 120051 : iGandhi et all l2006t ). The latest de- 
termination of w(9) from the XMM-Newton Large Scale 
Structure Survey gave u>(30") ~ 0.2 for ~ 1130 sources 
detected over 4.2 d eg 2 on the sky in the 0.5-2 keV band 

II Gandhi et al . 2006), although statistical significance is only 
~ 2(7. This result may be interpreted as the bias being ~ 3.7 
(or < 3.7 at the 2-cr level), which we have obtained as fol- 
lo ws: using the LDDE X LF and the best-fitting parameters 
in lHasinger etlH (|2005T ). we have calculated w(9) from equa- 
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tion Q, in which we have used the spatial correlation func- 
tion of non-linear dark matter fluctuations, as non-linearity 
cannot be ignored at such a small a ngular separ ation. (We 
have used the halo model approach (|Seliakll2OO0l ) to obtain 
the non- linear power spectrum.) If such a large bias is re- 
alized also for blazars, then GLAST should quite easily be 
able to detect the correlation. We should, however, keep in 
mind that no significant correlations were observed in the 
2-10 keV band, and one needs to improve the source statis- 
tics before making any definitive conclusion about the bias 
of AGNs from X-ray surveys. 



3.2.3 Bias from halo model 

Finally, we try to estimate the bias of blazars from the halo 
model. As any galaxies (hosting blazars) must form in dark 
matter haloes, the bias of blazars should be related to the 
bias of dark matter haloes, bh(M,z), which is known accu- 
rately from iV-body simulations as well as from analytical 
models such as th e extended Press-Schechter model (e.g., 
iMo fc Whitel fl996). Let us suppose that the 7-ray lumi- 
nosity of blazars is correlated with the mass of host dark 
matter haloes, L 7 = Lj(Mh). Then, the bias of blazars 
may be estimated from the bias of dark matter haloes via 
b B (L~,,z) = bh(Mh(L~ ,),z). We use b h (M h (L^), z) derived 
by IMo fc White] (| 19961 ). 

The relation between the 7-ray luminosity of blazars 
and the host halo mass is not known. Whether or not such a 
relation actually exists is not known either. Moreover, even if 
there is a relation between L 7 and Mh, not all haloes would 
host blazars, especially when one takes into account the fact 
that the jets from blazars should be directed towards us. We 
nevertheless estimate bs here using the following argument. 
Since it is plausible that the blazar 7-rays are emitted via the 
inverse-Compton scattering of the relativistic electrons ac- 
celerated in blazar jets off the cosmic microwave background 
(CMB), the 7-ray luminosity would be correla ted with the 
activit y of a central supermassive black hole. IWang et al.l 
(2002) have found an empirical relation between the blazar 
luminosity at the peak frequency (i.e., a frequency at which 
the blazar emits most energy, vL v ), L p k, in the optical to 
X-ray regime, and the luminosity of emission lines, Lii ncs , in 
the optical. The line luminosity is related to th e Eddington 
lumin osity via log(Lii ncs /I/Edd) ~ —5 to —3 l)Wang et al.l 
2002) which, in turn, would give a L p ^-Mbh relation, where 
Mbh is a black hole mass which determines the Eddington 
luminosity. We then relate the 7-ray luminosity, L 7 , with the 
peak luminosity, L pk, by assuming that they are roughly 
equa l, L 7 w L p k (|lnoue fc Takaharal 1 19961 ; iKataoka et al.l 
1 19991 ). We still need to relate the black hole mass, Mbh, 
with the host halo mass, Mh, for which we use a correla- 
tion betweeji_thel2lackhole mass and the host galaxy mass, 
M (|Ferrarese et al.l l2006). Of course, the galaxy mass may 
not be equal to the dark halo mass, the latter being larger, 
as dark matter haloes extend more than the luminous part 
of galaxies. We find that this uncertainty hardly affects our 
prediction for the blazar bias. Using these arguments, we 
finally obtain the desired relation, 

M = 1011 ^( ib^g^r ) 17 ' ™ 




Redshift z 

Figure 7. The same as Fig. [5] but for the blazar bias equal to 
the bias of dark matter halos with an empirical 7-ray luminosity- 
mass relation [equation (122 H . 



for log(I/ii nes /LBdd) = —4. We shall use this relation for 
estimating the blazar bias from 6s(L 7 ,z) = bh(M (L 7 ) , z) . 

Before we proceed further, let us address the uncer- 
tainty in our treatment. If we used log(Z/rm es /-I<Edd) = —3 
instead of —4, and L 7 = 0.1L p k instead of L 7 = -Lpk, we 
would obtain M = 1O 13 M (L 7 /1O 44 ergs" 1 ) 17 , making 
the host haloes more massive and hence the larger bias for 
the same luminosity. We find, however, that the resulting bs 
did not change significantly. 

Figure |7Ja) shows &b(«) from equation J5]) with 
b B {L 1 ,z) = b h (M(L 1 ),z) and equation ([22]). Figure JT^b) 
shows that the most dominant contribution to w(8) at 
8 = 20° comes from z < 0.1, which implies that the an- 
gular power spectrum of blazars that would be detected by 
GLAST would have the average bias of about 0.4, which 
is too small for GLAST to measure Cf . Both LDDE1 and 
LDDE2 give very similar results. The reason why we found 
such a small bias is that our model predicts that GLAST 
would detect fainter blazars at low redshifts than bright 
galaxies at high redshifts. Therefore, the average bias is dom- 
inated by the faint, low-z blazars that have a small bias. 
(Faint blazars are formed in low-mass dark matter halos, 
which have a small bias.) At 2 > 0.5, on the other hand, 
5b (z) can be as larger as unity, as only the bright blazars 
(in massive haloes) are detectable. 

Although it does not affect our analysis very much, let 
us mention one subtle feature of the halo approach we just 
described. The GLF of blazars may be related to the mass 
function of dark matter halos, dnh/dMh, as 

dL 7 Pl {L 1 ,z)=dM h ^-(M h ,z)N B {M h ,z), (23) 
dMh 

where Nb (Mh , z) is the so-called "halo occupation distri- 
bution," which represents the average number of blazars 
per each halo of mass, Mh, at a given redshift, z. One 
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can use this relation, GLF (p T ), and L^—Mh relation 
to obtain NsjMh, z). U sing the Press- Schechter function 
l|Press fc Schechteri 1 1974 ) for dn h /dMh and the LDDE1 
model for p 7 , we have found that Nb(M,z) diverges ex- 
ponentially at the high- mass end. This is because dn^/dMh 
has an exponential cutoff, while p 7 given by the LDDE1 
model decreases only as a power law with the luminosity. 
Of course this divergence is an artifact from the fact that 
we do not know the precise shape of the GLF at the bright- 
est end, which is poorly constrained. It is likely that (i) the 
GLF has a maximum luminosity above which it rapidly ap- 
proaches zero, and (ii) the assumed L 7 -M relation cannot 
be extrapolated to large luminosities. In order to remove the 
divergence, we have put an upper cutoff in the GLF so that 
Nb(M, z) never exceeds 1. For the LDDE1 model the cutoff 
luminosity is 3 x 10 47 erg s" 1 , which keeps the GLF still 
consistent with the EGRET data because the contribution 
to the GLF fr om such luminous b l azars is not significant 
(see Fig. 12 of iNarumoto fc Totanil i|2006r )i 

Any of the assumptions we have made in this subsection 
could be incorrect. The GLAST data will provide us with 
much better idea about the clustering of blazars, which will 
enable us to test these assumptions. 



4 ANGULAR POWER SPECTRUM OF 
GALAXY CLUSTERS 

4.1 Gamma-ray emission from galaxy clusters 

There is a fascinating possibili ty that GLAST finds 



clusters of galaxies in 
Colafrancesco fc Blasil Il99: 



7-ra ys ([Bcrczmsky et al.l 
lLoeb fc Waxmanl 
Waxman fc Loebl 



Totani & Kitavama 2000; _ 
Keshet et al.l 120031 . l2004al : iMiniatil 120031 ). A fraction of 



the EGRB may be due to these clusters of galaxies. 
While no convincing detectio n has been made so far from 
EGRET (|Reimer et ail 120031 ). there are a few reports of 
marginal evidence for correlatio n between the position of 
clusters and the EGRET data (|Kawasaki fc Totanil [20021 ; 
IScharf fc Mukheriedl20"02T ). In fact, there are many EGRET 
unidentified sources whose positions are coincident with 
Abell clusters or high galaxy density regions, but the 
physical association cannot be established because of low 
statistics and large EGRET error circles. It is expected that 
GLAST will give us the first conclusive evidence for 7-ray 
emission from clusters. 

Clusters of galaxies may emit 7-ray via two pro- 
cesses. One is the collision between relativistic protons ac- 
celerated by shock waves and surrounding cold matter, 
mainly protons, producing neutral pions which decay into 
7-rays. Since protons hardly lose their energy by radia- 
tive loss and their diffusion time is much longer than the 
age of the universe, all the cluster s are expected to emit 
some 7-rays by this mech anism l|Berezinskv et ail 1 19971 : 
IColafrancesco fc Blasil [l998h . The other mechanism is the 
inverse-Compt on scattering of relativi s tic electrons off the 
CMB photons jLoeb fc Waxmanll2000l; iTotani fc Kitavamal 



(JMtJ photons llboeb fe W axman 2UUU; lot am & Kitavama 
20001; IWaxman fc Loebl |2000| ; iKeshet et ail 120031 . 12004a 



Miniatill2003T r As for the source of these relativistic elec- 



electrons to relativistic speed. In the following sections we 
shall explore these two possibilities and calculate the angular 
power spectrum that would be measured by GLAST. 



4-1.1 Proton-proton collisions 

The cumulative luminosity function, <f>(z) [equation JBJ], for 
the cluster 7-ray emission from proton-proton collisions is 
given by 



f>C, PP (z) = 



dM h 



M h (F 7 , Um ,z) 



dn h 
dM h 



(M h ,z), 



(24) 



trons, the most popular scenario is that the shocks associ- 
ated with the formation of large-scale structure accelerate 



where we label quantities regarding clusters by attaching the 
subscript C henceforth; another subscript pp means that the 
7-ray emission comes originally from the proton-proton (pp) 
collisions. 

In order to relate the 7-ray flux to the halo 
mass and redshift, MhjF^^ z), we follow the model of 
IColafrancesco fc Blasil (|l998i l: relativistic protons are in- 
jected from the very central region of clusters, in which 
the central AGN or cD galaxy powers such an injection. 
These protons then diffuse from the central region to out- 
side with the efficiency that is determined by the magnetic 
field strength. We assume that a fraction, es = 10 -3 , of 
the baryon energy is given to the magnetic energy. For pro- 
ton energies of our interest, the diffusion time scale is always 
longer than the age of the universe; thus, protons are always 
confined within clusters. Using the radial injection profile of 
these diffused protons as well as the density profile of the 
surrounding medium that is well measured in X-rays, one 
can compute the rate of pp collisions. The effi ciency of the 
7-ray production from each collision is given in lKelner et ail 
(2006), which we follow in our calculation. We use a power 
law with an index of a p = 2.2 for the proton spectrum, 
with an upper cutoff whose energy is determined by a bal- 
ance between the diffusion time scale and the cluster age 
l|Colafrancesco fc Blasi|[l99^ ) . (We found that the cutoff en- 
ergy is much larger than the energy scale of our interest.) 

We calculate the total energy of relativistic protons, E p , 
by assuming that a fraction, e p , of the gravitational binding 
energy of baryons is given to protons, i.e., E v — e p Et ~ 
e p (Q.b/il m )MVc , where V c is the circular velocity at the 
virial radius, and adopt three values for e p = 0.5 (ppl), 0.1 
(pp2), and 0.01 (pp3). Note that the equipartition model, 
ppl, has been excluded marginally by observations (|Blasil 
ll999T ); however, we keep this model as an upper bound on the 
7-ray emission from clusters of galaxies via proton-proton 
collisions. The other two models (pp2 and pp3) are allowed 
by observations. The most pessimistic case s, pp3 and IC3, 
are in agreement with the estimates given in lGabici fc Blasil 
(2004). We summarize the parameters of these models as 
well as the expected number of clusters that would be de- 
tected by GLAST in Tabled We find that a large number of 
galaxy clusters are expected to be seen in the GLAST data, 
which would provide an exciting possibility of investigat- 
ing the physics of galaxy clusters using 7-ray observations. 
Figure [8ja) shows the redshift distribution of clusters that 
would be detected by GLAST with these models. 
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Table 2. Model parameters, a and e, the expected number count, 
N, the surface density, jV, the average bias, be (z = 0.01), of 
clusters that would be detected by GLAST. The last column lists 
the expected signal-to-noise ratio for detecting the correlation 
power spectrum averaged over 2 ^ I sC 30. 



Model 



/V 



AT (sr- 1 ) b c 



c°/sc 



ppl 


2.2 
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Figure 8. Redshift distribution of clusters of galaxies that would 
be detected by GLAST for (a) proton-proton collision and (b) 
inverse-Compton scattering models. Model parameters are given 
in Table [2] 



4-1.2 Inverse-Compton scattering 

Since electrons lose their kinetic energy via radiation 
loss rapidly compared with the dynamical time of clus- 
ters of galaxies, 7-ray emission would emerge only near 
the formation of shocks. The cumulative luminosity func- 
tion, 4>{z) [equation iJBJ], for the cluster 7-ray emission 
from the inverse-Compto n scattering is thus given by 
l|Totani fc Kitavamallioool ) 



0c,ic(z) 



L 



dM h R toxm (M h ,z)At 1 



(25) 



M h (F y>Um ,x) 



where Rf OIU1 (Mh, z) is the formation rate of clusters with 
mass of Mh at z, per comoving volume, and At 7 is the time 
scale during which 7-rays are radiated efficiently from each 
cluster. We calculate At 1 as either the inverse-Compton 
cooling time or the shock wave propagation time (whichever 
is longer): Ai 7 = max{tic, tshock}- In most cases of our in- 
terest, the latter is always much longer than the former, and 
therefore, Ai 7 = i s hock — r v i r /v 3 = 1.5(1 + z)~ 3 ^ 2 Gyr, in- 



dependent of Mh, where r v i r is the virial radius and v s the 
sound speed. The formation rate of clusters, R{ OTU1 (Mh, z), 
is given by the time-derivative of the halo mass function, 
dnh fdMh(Mh , z), correc ted for the halo destruction rate 
ijKitavama fc Sutol[l99r3 ). 

Similar to the proton-proton collision case, we calcu- 
late the total energy of relativistic electrons, E e , by as- 
suming that a fraction, e e , of the gravitational binding 
energy of baryons is given to electrons. We use a power 
law with an index of a e (either 2 or 2.2; see Table [2} for 
the 7-ray spectrum, with an upper cutoff whose energy 
is determined by a balance between the acceleration time 
scale and the cooling time scale. To calculate the acceler- 
ation time scale we use the magnetic field energy given by 
6b = 10 -3 times the binding energy of baryons. We choose 
(a e ,e e ) = (2,0.05), (2,0.01), and (2.2,0.01) as our models, 
and we call them IC1, I C2, and IC3, respect i vely. T he IC1 
model is investigated bv lTotani fc Kitavamal l|2000j) . and it 
gives maximally allowed number of 7-ray emitting clusters, 
as the IC1 model predicts the EGRB flux that is as large 
as what is measured by EGRET. These models are again 
summarized in Table [2] and Fig. [8] 



4.2 Angular power spectrum of galaxy clusters 
from GLAST 

The angular power spectrum of clusters of galaxiefl is given 
by equations <[U — (SJ) with the averaged blazar bias, bs, re- 
placed by the average cluster bias, 



bc, PP (z) = 



d M h ^-(M h ,z) 



<t>C,pp(z) J A f h (F 7 lim ,z) 

xb h (M h ,z), (26) 

1 f°° 

&c,ic(z) = 7— r / dM h R{ rm{M h ,z)At-, 

(27) 



x b h {M h ,z), 



for the proton-proton collision model and the inverse- 
Compton model, respectively. 

Figure [9] shows the angular power spectrum of these 
7-ray clusters with the binned error boxes (Al — 0.5Z) as 
well as the shot noise term for (a)-(c) proton-proton colli- 
sion and (d)-(f) inverse-Compton models. The correlation is 
quite significant particularly for optimistic models predict- 
ing large number of 7-ray emitting clusters being detected 
by GLAST, i.e., ppl and IC1. The last column of Table [2] 
shows the signal-to-noise ratio for the power spectrum av- 
eraged over 2 < I < 30, C/SC. We find that the signal-to- 
noise ratio exceeds unity for all the models that we have 
considered: the minimum is C/SC = 1.3 for IC3, and the 
maximum is 7.0 for ppl, despite the fact that only small 
number of clusters are expected to be seen in the GLAST 
data. This is because clusters of galaxies are formed in the 
high-density peaks and thus are highly biased. The sixth 
column of Table [2] shows the average bias factors of clusters 
at z = 0.01. 



4 IWaxman fc Loebl ||2000|) , followed by iKeshet et all ||2003| . 
l2004al ). studied the angular correlation of the radio and 7-ray 
background radiation from galaxy clusters. 
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Figure 9. The same as Fig. [3]but for galaxy clusters. The left panels show the proton-proton collision models (ppl, pp2, pp3), while the 
right panels show the inverse-Compton models (IC1, IC2, IC3). See Table[2]for the model parameters. 



5 DISCUSSION 

5.1 Admixture of blazars and galaxy clusters 

While follow-up programs should reveal the identity of the 
GLAST 7-ray sources and also some of the galaxy clus- 
ters might appear as extended sources, at very early stage 
of GLAST observational campaign, all the point sources 
should more generally be considered to be mixed of various 
emitters. Here we consider two-population case, blazars and 
galaxy clusters. Our purpose in this section is to investigate 
whether it is possible to distinguish the blazar component 
from that of clusters by the angular clustering, even before 
the follow-ups. 

When there are more than one species of sources on 
the sky, one should also consider cross-correlation between 
different species. When there are blazars and galaxy clusters 
in the 7-ray sky, the angular power spectrum is given by 



Ci = C;,b + Ci,c + ZC^bc, 



(28) 



where C;,b and C;,c are the spectra from blazars (Section[2]) 
and clusters (Section [4}, respectively. The surface number 
density of sources in this case is instead given by the sum 
of the two species, N = Mb + Mc- The last term, Ci,bc, 
represents the cross-correlation between blazars and galaxy 



clusters, and is given by 



C,, BC = 2tt / dcosfl P l {cose)w BC {0), (29) 



N 2 w BC {0) = 



dz ^L X{Z)2 



X(t>B{z)4>c{z)b B {z)bc(z) 
> I du ^(V^ + xO^v) , (30) 



where <^>b(z) [equation ©J and 4>c{z) [equation (|24|l or (I25[) ] 
as well as 1>b{z) [equation ([5]l] and be (z) [equation (|26[) or 
(12711 ] have been given in the previous sections. 

Figure [10] shows the angular power spectrum of GLAST 
point sources, including blazars (LDDE1 and 6s = 3) and 
galaxy clusters for various models. If the average blazar bias 
is as large as 3, the power spectrum from blazars almost 
always dominates the signal, especially at low Vs. An inter- 
esting feature in Fig. [TJJ is that the shape of blazar power 
spectrum and cluster spectrum are quite distinct: the clus- 
ter spectrum falling towards low multipoles more rapidly 
than the blazar spectrum (i.e., the cluster spectrum rising 
towards high multipoles more rapidly than the blazar spec- 
trum). Therefore, when the blazar signal dominates (pp3 and 
IC3 in Fig. I10[) . one would see a shallower power spectrum, 
while when the cluster signal dominates (ppl and IC1), one 
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Figure 10. Angular power spectrum of blazars for the LDDE1 model with 6g = 3 (dotted), clusters of galaxies (dashed), and cross- 
correlation between blazars and clusters (dot-dashed). The thick solid curves show the total signal, Cf = Cf B + Cp c + 2C;bc, while 
the boxes show the binned errors (A( = 0.5i). The left panels show the proton-proton collision models (ppl, pp2, pp3), while the right 
panels show the inverse-Compton models (IC1, IC2, IC3). See Table[2]for the model parameters. 



would see a steeper power spectrum at / < 30. This feature 
may help us identify the dominant source of clustering seen 
by GLAST. We do not know which point sources GLAST 
would detect, blazars or clusters, until the follow-up observa- 
tions are carried out; however, the angular power spectrum 
may provide us with useful information regarding the dom- 
inant species. 

5.2 Identifying GLAST point sources with radio 
survey 

The Faint Images of the Radio Sky at Twenty Centimeters 
(FIRST) survey (jBecker et all Il995l ) has detected ~ 811k 
radio sources over 9,033 square degrees, or ~ 90 sources per 
square degree, with the source detection threshold of 1 mjy 
at 1.4 GHz0 Since many of the GLAST point sources will 
be within the FIRST field of view, the FIRST survey will 
provide us with valuable information regarding identification 
of the point sources that would be detected by GLAST. 

How bright are blazars and galaxy clusters in radio? 
The EGRET blazar catalogue was constructed such that 

5 http://sundog.stsci.edu/first 



the EGRET sources are also detected in radio. The cor- 
relation between radio and 7-ray luminosities of blazars is 
not yet well understood , and there is a considerable disper- 
sion |Miicke et al.ll 19971 ). However, the standard synchrotron 
self-Compton model of blazars predicts t hat t here should 
be some correlation. iNarumoto fc Totanil |2006l) assumed a 
proportional relation between mean radio luminosity, L r , 
and L 7 , with a dispersion obeying to the log-normal dis- 
tribution, and found that mean relation of L r = 10 _3 ' 23 L 7 
and dispersion er(rog(L r / X 7 )) = 0.49 fitted the observed 
data well. Using this relation and assuming that the spec- 
tral index in radio is given by a r — 1.0 (because the radio 
emission is due to synchrotron radiation), we find that the 
flux in radio that corresponds to the limiting flux for the 
point sources that can be detected by GLAST in 7-rays is 
given by F P ,GLAST,ito = (1 + z) 2 ~ a - (10- 3 ' 23 L 7 )/(47rd 2 ) = 
10~ 323 (1 + z) a ''~ ar (a 1 - l)-E min F 7i i im , where we have used 
equation l|15[l. We thus find i*V,GLAST,iim ~ 10 mjy, which is 
an order of magnitude brighter than the limiting flux for the 
point sources detected by the FIRST surveyQ Therefore, we 



6 The luminosities at 1.4 GHz and at 2.7 GHz are the same, as 
we adopt a r = 1. 
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expect the radio counterparts for the GLAST blazars to be 
found in the existing FIRST point source catalogue (if they 
are within the FIRST field of view), although some sources 
which deviate from the L r -L 1 relation above may be missed. 

The radio emission from galaxy clusters is also likely 
from synchrotron radiation. For the inverse-Compton sce- 
nario one can estimate the luminosity in radio from a ratio 
of the CMB and magnetic field energy density, in which case 
the radio luminosity from galaxy clusters is much smaller 
than that expected from blazars ha ving the same 7-ray lu- 
minosity (jTotani fc Kitavamall200oh . For the proton-proton 
collision scenario they would be even dimmer in radio; other- 
wise they should also be detectable by 7-rays due to the elec- 
tron inverse-Compton scattering. Therefore, unlike blazars, 
galaxy clusters would not be identifiable with the FIRST 
survey, which makes the FIRST survey a good diagnosis 
tool for identification of the GLAST point sources. 

5.3 Measuring blazar anisotropy with GLAST 

Based upon the results that we have obtained so far, we 
here show one example strategy for the point source survey 
and identification of blazars that would be carried out by 
GLAST. 

(i) Source detection. — After its launch, GLAST will 
start detecting 7-ray sources from all the directions on the 
sky. Some of them would be extended sources (such as 
nearby galaxy clusters), and some would be highly variable 
(such as 7-ray bursts). These sources should be removed. 

(ii) Removing galaxy clusters. — As blazars are also 
bright in radio but galaxy clusters are not, one may remove 
galaxy clusters from the GLAST data using the source cat- 
alogue from the FIRST survey (Section l5.2p . 

(iii) Updating GLF and further cut. —With the GLAST 
source catalogue from (ii), which would consist mostly of 
blazars, one may update the GLF of blazars by extending 
it down to fainter sources than those that EGRET has de- 
tected. At this point we probably gain some insight as to 
which GLF model fits the data better, LDDE1 or LDDE2, 
or whether or not we need a different GLF model. If LDDE1 
is indeed confirmed, then one needs other populations of 
sources in order to explain the bulk of EGRB. 

(iv) Analysis of angular power spectrum. — Measure the 
angular power spectrum of the sources that have survived 
the cuts in (i) and (ii). Since we have currently several mod- 
els which predict a variety of the blazar bias, from 0.4 to 4 
(Section l3.2|l . the power spectrum measured at I < 30 should 
provide us with useful information regarding the blazar bias. 
While we would expect the contribution from galaxy clusters 
is minimal in this catalogue owing to the cuts in (i) and (ii), 
the shape of the power spectrum would also provide useful 
(albeit indirect) confirmation that the bulk of the sources 
in the catalogue are blazars. The blazar bias measured from 
GLAST, or an upper limit on the bias, would be the first 
direct measurement of the bias of blazars, which would shed 
light on the formation process of blazars and their link to 
the quasars detected in the optical and the AGNs detected 
in X-rays. 

(v) Completion of follow-ups: beginning of precision 
study. — When the source identification with direct follow- 
up observations is complete, one should revisit the blazar 
source catalogue again, establish the GLF of blazars more 



firmly, and re-analyze the angular power spectrum. At this 
point it would also be possible to obtain the 3-d power spec- 
trum, as opposed to the angular spectrum, using the red- 
shift information from follow-up observations. This would 
be very powerful in constraining the formation and evolu- 
tion of blazars, as one can constrain the evolution of blazar 
bias as a function of redshift, provided that enough number 
of blazars are detected by GLAST. 



6 CONCLUSIONS 

In this paper we have calculated the angular power spec- 
trum of blazars and galaxy clusters that would be detected 
by GLAST. We have shown that GLAST can detect the 
spatial clustering of blazars if the average bias of blazars 
exceeds 1.2 and 0.5, for the canonical GLF model (LDDE1) 
that accounts for 25-50% of the extragalactic 7-ray back- 
ground (EGRB) and the extreme mod el (LDDE2) that ac- 
coun ts for all the EGRB, respectively l|Narumoto fc Totanil 
2006). While the blazar bias is not known with any preci- 
sion, current observations seem to suggest, albeit indirectly, 
that it can take on any values between ~ 0.4 and ~ 4; thus, 
the GLAST data will provide us with the first, direct esti- 
mate of the bias of blazars which, in turn, would constrain 
the formation and evolution of blazars. 

As for galaxy clusters, which are highly biased ob- 
jects, we have found that the signal-to-noise ratio of the 
correlation exceeds unity for all the models we have con- 
sidered (TalDky2]h_pj^5tor^^ 

decay (Berezinsky et al 1 1 19971 ; IColafrancesco fc Blaslll998l ). 
and in verse-Compton scatterin g of relativistic electrons off 
CMB (lLoeb fc Waxmanl |2000|; iTotani fc Kitavamal l200d; 
IWaxman fc Loebl |2000| ; iKeshet et all 120031 . l2004al ; iMiniatl 
l2003h . 

We have shown that the angular power spectra of 
blazars and galaxy clusters are quite distinct at low mul- 
tipoles, I < 30, the blazar spectrum being much shallower 
than the cluster one. This feature helps us identify the pop- 
ulation dominating the angular power spectrum of the point 
sources that would be detected by GLAST. 

Although the full follow-up observations would take 
long time, a quick (but less accurate) identification of sources 
is possible with the existing FIRST survey data in radio at 
1.4 GHz, as most of the blazars should also be bright enough 
to be seen in radio, while galaxy clusters should not. We 
have given an example strategy for using the angular power 
spectrum as a diagnosis tool for blazars in Section \5. 3 1 

With an impressive number of blazars as well as galaxy 
clusters expected to be detected by GLAST, one should 
maximize the scientific outcome from the GLAST data by 
using as many tools as possible. The angular power spec- 
trum (or angular correlation function) is easy to calculate 
from the point source catalogue, and would give us invalu- 
able information about the spatial clustering of blazars and 
high-energy activity in clusters of galaxies, which are poorly 
known at present. 
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